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Abstract.  For  a three-layer  system  consisting  of  two  ferromagnetic  layers  separated  by  an  antifer- 
romagnetic interlayer  it  is  shown  that  the  stability  region  of  single-domain  ferromagnetic  layers  is 
technologically  dependent  on  the  ratio  between  the  width  of  interface  atomic  steps  and  the  thickness 
of  layers.  The  thickness-roughness  phase  diagram  has  three  regions:  a collinear  orientation  of  the 
magnetizations  of  ferromagnetic  layers  (1),  oriented  to  each  other  at  90°  (2),  and  a multidomain 
structure  of  ferromagnetic  layers  (3).  The  proposed  approach  can  be  compared  with  experiment  by 
studying  magnetic  hysteresis  loops  against  nanostructure  layer  thicknesses  and  their  roughnesses 
owing  to  atomic  steps  on  interfaces. 


Introduction 

The  discovery  of  giant  magnetoresistance  [ 1 ] has  stimulated  interest  in  sandwiches  consist- 
ing of  alternating  thin  ferromagnetic  (Fe,  Co)  and  nonmagnetic  (Cr,  Cu)  metallic  layers. 

An  antiparallel  or  parallel  orientation  of  adjacent  ferromagnetic  layer  magnetization  is 
energetically  favorable  when  the  number  of  atomic  planes  in  the  antiferromagnetic  spacer 
is  even  or  odd,  as  in  Fig.  1(a)  and  1(b)  respectively.  In  the  latter  case  an  external  magnetic 
field  changes  the  antiparallel  orientation  of  magnetization  to  parallel  one.  It  is  accompained 
by  a drop  in  the  electric  resistance  up  to  tens  of  percent. 

An  example  of  such  a layered  antiferromagnet  is  a chromium  film  with  the  thickness 
32  A < d < 150  A,  in  which  at  low  temperature  there  appears  a commensurate  transverse 
spin-density  wave  [2],  A similar  structure  occurs  when  the  iron  impurities  with  the  con- 
centration exceeding  2%  are  introduced  into  chromium  [3], 

The  roughness  of  the  interfaces  owing  to  atomic  steps  on  them  can  make  the  uniform 
order  parameter  distribution  energetically  unfavourable,  because  the  spin  orientations  of 
adjacent  atoms  on  the  step  are  opposite  (Fig.  2).  If  the  characteristic  distance  R between 
atomic  steps  exceeds  some  critical  value,  it  is  energetically  favourable  to  break  up  the  layers 
into  domains  with  parallel  and  antiparallel  orientation  of  magnetization  [4]. 

In  this  paper  we  report  a phase  diagram  of  multilayer  ferromagnet-antiferromagnet 
nanostructures  with  arbitrary  ratios  between  the  exchange  interactions  as  well  as  between 
the  layer  thickness  and  the  characteristic  separation  of  steps  at  interfaces.  The  proposed 
approach  can  be  used  in  the  experiments  studying  the  forms  of  magnetic  hysteresis  loops. 

1.  Description  of  the  model 

We  shall  investigate  the  disrtibution  of  order  parameters  in  layers  in  the  mean-field  approxi- 
mation [5,6].  The  parameters  can  be  introduced  for  each  layer  as  follows:  the  magnetization 
vector  for  ferromagnetic  layers  and  the  antiferromagnetism  vector,  equal  the  difference  in 
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Fig.  1.  Orientation  of  spins  in  a three-layer 
system  consisting  of  two  ferromagnetic  layers 
and  an  antiferromagnetic  interlayer  in  the  case 
of  smooth  interfaces  and  odd  (a)  and  even  (b) 
number  of  atomic  planes  in  the  interlayer. 


Fig.  2.  Orientation  of  spins  near  the  step  at 
the  interface  for  homogeneous  distribution  of 
the  order  parameters  (a)  and  in  the  presence  of 
a domain  wall  (b). 


the  magnetizations  of  sublattices,  for  antiferromagnetic  layers.  In  case  of  sufficiently  thin 
layers,  the  spins  of  atoms  lie  in  the  plane  and,  therefore,  the  orientation  of  the  vector  order 
parameter  can  be  given  by  the  angle  8 that  it  forms  with  the  x axis,  which  lies  in  the  plane 
of  the  layer. 

Far  from  Curie  and  Neel  temperature,  the  exchange  energy  W,-  caused  by  the  interaction 
inside  the  ;th  layer  can  be  written  as 


w<  = Y f (V0i)2  dv’  (1) 

where  the  integral  is  taken  over  the  volume  of  the  layer,  and  A,-  is  the  corresponding 
exchange  constant.  In  the  order  of  magnitude.  A,-  ~ ./,  Sj / b , where  7;  is  the  exchange 
integral  between  adjacent  atoms,  b is  the  interatomic  distance,  and  5,-  is  a mean  spin  of 
an  atom. 

The  energy  of  the  exchange  interaction  between  adjacent  layers  with  numbers  i and 
i + 1 can  be  written 

WiJ+ 1 = ±B  J cos(0j  -ft+i)dS,  (2) 

where  integration  is  performed  over  the  layer  interface,  B ~ J f,afSiSi+i/b2 , J f,af  is  the 
exchange  integral  between  adjacent  atoms  belonging  to  different  layers,  and  the  sign  " + " 
is  opposite  for  different  sides  of  an  atomic  step  at  the  interface  of  the  layers. 

We  divide  all  distances  by  b,  assuming  it  to  be  virtually  identical  in  both  sorts  of  layers 
and  divide  all  energies  by  the  constant  Aaf  for  an  antiferromagnetic  layer.  We  introduce 
dimensionless  parameters 


a — 


Jf,afSf  a JfS2f 

, ' and  y — 

JafSaf  JafSaf 


(3) 


where  the  subscripts  / and  af  correspond  to  a ferro-  and  antiferromagnet. 

By  varying  the  parameters  0/  in  (1),  we  obtain  a differential  equation  that  describes  the 
distribution  of  the  order  parameter  in  a layer 


Mi  = 0. 


(4) 
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A more  careful  procedure  is  needed  to  obtain  correct  conditions  at  the  interface  of  layers. 
The  energies  of  intra-  and  interlayer  interactions  must  be  varied  in  a discrete  model,  and 
then  we  must  go  to  a continuum  limit.  As  a result  we  get 

ddi  B 

A 0;  - — = Ttt  sin  ((9,-  - <9/+1),  (5) 

dn  Ajb 

where  A is  a two-dimensional  Laplacian  in  the  plane  of  a layer,  d/dn  is  the  derivative  in 
the  direction  of  an  o uter  normal  to  the  boundary  of  the  layer,  and  upper  and  lower  signs  in 
(5)  correspond  to  those  in  (2). 

Thus,  to  find  the  distribution  of  order  parameters  in  a multilayer  nanostructure,  it  is 
necessary  to  solve  the  system  of  linear  differential  equations  (4)  with  non-linear  boundary 
conditions  (5).  The  distribution  will  depend  on  the  values  of  a and  y,  on  the  thickness 
of  layers,  as  well  as  on  the  characteristic  distance  R between  the  steps  at  the  interface  of 
layers. 

2.  Phase  diagram  of  a three-layer  system 

The  results  of  model  calculations  can  be  shortly  displayed  on  the  phase  diagram  (Fig.  3). 
Phase  1 is  characterized  by  the  presence  of  static  vortices  at  the  interfaces  and  a collinear 
orientation  of  magnetization  of  ferromagnetic  layers.  In  phase  2,  the  magnetizations  of 
the  ferromagnetic  layers  are  homogeneous  and,  in  the  absence  of  an  external  magnetic 
field,  are  oriented  to  each  other  at  90°.  Phase  3 corresponds  to  a multidomain  structure  of 
ferromagnetic  layers. 

The  phase  diagram  can  be  compared  with  experimental  data  by  studying  the  state  of 
ferromagnetic  layers  using  a magnetic-force  microscope  with  various  ratios  between  R and 
a and  different  temperatures,  because  a and  y are  temperature-dependent. 


3.  Effect  of  a magnetic  field:  hysteresis  loops 


The  magnetic  flux  reversal  of  ferromagnetic  layers  in  phase  1 occurs  virtually  independently, 
and  the  hysteresis  loops  must  coincide  with  those  in  a two-layer  system  consisting  of 
one  ferromagnetic  and  one  antiferromagnetic  layer.  Here  and  below  we  assume  that  the 
maximum  magnetic  field  is  much  less  than  the  exchange  field  in  the  antiferromagnet. 
Therefore,  the  magnetization  of  antiferromagnetic  layers  can  be  neglected. 


Fig.  3.  Thickness-roughness  phase  diagram 
for  a three-layer  system  consisting  of  layers 
of  identical  thickness.  The  region  of  weak 
distortion  of  order  parameters  is  hatched. 


Fig.  4.  Hysteresis  loops  of  a three-layer  nanostruc- 
ture consisting  of  two  ferromagnetic  layers  sepa- 
rated by  an  antiferromagnetic  interlayer  with  iden- 
tical thickness. 
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In  phase  2,  in  a weak  magnetic  field  that  exceeds  the  anisotropy  field  in  the  plane  of 
ferromagnetic  layers,  the  magnetizations  of  ferromagnetic  layers  are  oriented  at  45°  to  the 
field  while  remaining  virtually  perpendicular  to  each  other.  In  this  case,  the  magnetization 
of  the  system  is  Mmax/-s/ 2,  where  Mmax  is  the  maximum  magnetization  of  ferromagnetic 
layers.  Its  further  evolution  can  be  studied  by  minimizing  the  sum  of  the  energies  of 
ferromagnetic  layer  interaction  with  each  other  and  with  the  magnetic  field.  The  energy  of 
a ferromagnetic  layer  in  the  magnetic  field  with  inductance  B is  given  by 

W f = — 2,/j.l  B cos  — , (6) 

where  /<  is  the  magnetic  moment  of  an  atom  of  a ferromagnet.  For  the  angle  4*  between 
the  magnetization  of  layers  we  have  a transcendental  equation 

= HlB  cos  j.  (7) 


The  characteristic  field  B*  (Fig.  4 ) at  which  there  occurs  a substantial  magnetization  change 
is 


ixal 


(8) 


In  phase  3,  in  a weak  magnetic  field,  the  domains  with  parallel  orientation  of  the  mag- 
netization of  layers  are  oriented  along  the  field.  In  this  case  the  system  magnetization  is 
Mmax/  2-  The  magnetizations  of  ferromagnetic  layers  in  domains  with  antiparallel  orienta- 
tion in  a zero  field  behave  similarly  to  the  magnetizations  of  the  sublattices  in  the  volume 
of  an  antiferromagnet.  They  are  oriented  virtually  perpendicular  to  the  external  field. 

As  B increases,  the  angle  T*  between  them  decreases.  The  characteristic  value  of 
induction  B*  of  the  external  magnetic  field  in  which  the  angle  varies  substantially,  is 
found  by  the  method  similar  to  that  for  phase  2 and  is  given  in  order  of  magnitude  by  (8). 
Therefore,  the  form  of  the  hysteresis  loops  in  phases  2 and  3 differs  only  in  the  magnit  ude 
of  magnetization  in  weak  fields. 
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